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Aims. We present evolutionary calculations and colors for massive white dwarfs with oxygen-neon cores for masses 
between 1.06 and 1.28 Mq. The evolutionary stages computed cover the luminosity range from \og{L/LQ) ~ 0.5 down 
to -5.2. 

Methods. Our cooling sequences are based on evolutionary calculations that take into account the chemical composition 
expected from massive white dwarf progenitors that burned carbon in partially degenerate conditions. The use of 
detailed non-gray model atmospheres provides us with accurate outer boundary conditions for our evolving models at 
low effective temperatures. 

Results. We examine the cooling age, colors and magnitudes of our sequences. We find that massive white dwarfs are 
characterized by very short ages to such an extent that they reach the turn-off in their colors and become blue at ages 
well below 10 Gyr. Extensive tabulations for massive white dwarfs, accessible from our web site, are also presented. 
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1. Introduction 

White dwarf stars constitute the most common end- 
product of stellar evolution. The present population of 
white dwarfs thus contains precise information about the 
star formation rate of our Galaxy, as well as about its age, 
which is information that can be accessible from their mass 
and luminosity distributions, as long as evolutionary mod- 
els for the progenitor of white dwarfs and for the white 
dwarfs themselves are available. However, most of the in- 
formation about the birth and the evolution of the galactic 
disk is concentrated at the faint end of the white dwarf lu- 
minosity function, which is dominated by the contribution 
of massive white dwarfs (Dfaz-Pinto et al. 1994). It is also 
worth mentioning at this point that the MACHO collabora- 
tion in their first season reported a microlensing event with 
a duration of 110 days towards the galactic bulge (Alcock 
et al. 1995). For this particular event a parallax could be 
obtained from the shape of the light curve, from which a 
mass of l.Sto^e-^O then derived, indicating that the 
gravitational lens could possibly be a massive oxygen-neon 
(ONe) white dwarf or a neutron star. 

The interest in very low-luminosity white dwarfs has 
also increased since the MACHO team proposed that the 
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microlensing events towards the LMC could be due to a 
population of faint white dwarfs — see, however, Isern et 
al. (1998a), Torres et al. (2002), and Garci'a-Berro et al. 
(2004). Since ONe white dwarfs cool faster than the bulk 
of carbon-oxygen (CO) white dwarfs, it is reasonable to ex- 
pect that perhaps some of the events could be due to these 
elusive massive white dwarfs. Moreover, studies about the 
distribution of masses of the white dwarf population (Finley 
et al. 1997; Liebert et al. 2005) show the existence of a nar- 
row sharp peak near 0.6 Mq, with a tail extending towards 
higher masses, with several white dwarfs with spectroscop- 
ically determined masses within the interval comprised be- 
tween 1.0 and 1.2 Mq. 

On the other hand, theoretical evidence suggests that 
high-mass white dwarfs should have cores composed mainly 
of oxygen and neon — at least for non-rotating stars 
(Dominguez et al. 1996) — in contrast to average-mass 
white dwarfs, for which carbon-oxygen cores are expected. 
The existence of such massive white dwarfs has been sug- 
gested as the result of either binary evolution (Marsh et 
al. 1997) — see also Guerrero et al. (2004) — or of the 
evolution of heavy-weight intermediate-mass single stars 
(Ritossa et al. 1996; Garcia-Berro et al. 1997b; Iben et 
al. 1997, Ritossa et al. 1999). In particular, Garcfa-Berro 
et al. (1997b) found that, when the core mass of a 9 Mq 
white dwarf progenitor exceeds ~ 1.05 A^q, carbon is ig- 
nited off-center in semidegenerate conditions before reach- 
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ing the thermally pulsing phase at the AGB tip. As a result 
of repeated carbon-burning shell flashes that ultimately 
gave rise to carbon exhaustion, these authors found that 
at the end of carbon burning the star was left with an 
oxygen-neon core almost devoid of carbon. After consid- 
erable mass-loss episodes, the progenitor remnant is ex- 
pected to evolve into the central star of a planetary nebula 
and ultimately into a white dwarf with an oxygen-neon 
core. A possible observational counterpart of these ultra- 
massive white dwarfs would be the single massive white 
dwarf LHS 4033, which has a mass of ~ 1.32 Mq (Dahn 
et al. 2004). Other possible massive white dwarfs with 
oxygen-neon cores would be the magnetic white dwarf 
PG 1658+441 (Schmidt et al. 1992; Dupuis et al. 2003) 
— with a mass of ~ 1.31 Mq — the highly magnetic white 
dwarf RE J0317-853 (Ferrario et al. 1997), which has a 
mass of ^ 1.35 Mq, and the ultramassive white dwarf 
GD 50 (Dobbie ot al. 2006). 

The mass radius relation and the pulsational properties 
of massive white dwarfs have been the subject of recent 
theoretical work: Althaus et al. (2005) and Corsico et al. 
(2004), respectively. However, the evolution (namely, the 
cooling ages and colors) of massive white dwarfs consider- 
ing the core composition as predicted by the evolution of 
massive progenitor stars has not yet been studied in detail, 
in sharp contrast with the situation for standard CO white 
dwarfs for which very accurate cooling sequences do exist; 
see, for instance, Salaris et al. (2000) and references therein. 
A first attempt to describe the cooling of ONe white dwarfs 
was performed by Garcia Berro et al. (1997a) using a sim- 
plified cooling code, but, although the equation of state 
employed in this work was very detailed, the evolutionary 
calculations were rather simplistic and the adopted chem- 
ical profiles were a flat ONe mixture, without taking into 
accoimt the CO buffer on top of the ONe core that full 
evolutionary calculations predict. This paper is aimed at 
precisely filling this gap by presenting new evolutionary 
calculations for massive white dwarfs with oxygen-neon 
cores down to very low surface luminosities and effective 
temperatures. In addition, we present colors and magni- 
tudes for these stars on the basis of non gray model atmo- 
spheres. Detailed model atmospheres also provide us with 
accurate outer boundary conditions for our evolving mod- 
els. Our paper is organized as follows. In §2 we present our 
input physics. In §3 we discuss our evolutionary sequences. 
Finally, in the last section we summarize our findings and 
draw our conclusions. 

2. Input physics 

The evolutionary code and the starting white dwarf con- 
figurations used in this work are essentially those used re- 
cently in the calculation of mass-radius relations for mas- 
sive white dwarfs by Althaus et al. (2005), and we refer the 
reader to that paper for further details. Because we are now 
interested in providing accurate cooling times, some major 
improvements to the input physics assumed in Althaus et 
al. (2005) have been made. First, we have included the re- 
lease of latent heat upon crystallization. Despite the fact 
that crystallization in massive white dwarfs takes place at 
higher stellar luminosities than do standard carbon-oxygen 
white dwarfs, its impact on the cooling times is not entirely 
negligible. In our calculations, crystallization sets in when 
the ion coupling constant T= 180, where F = Z^e^/fkBT 



and r is the radius of the Wigner-Seitz sphere. We con- 
sidered a latent heat release of per ion, which was 
spread over the range 175 < F < 185. Phase separation 
upon crystallization has been shown to introduce negligi- 
ble time delays for massive white dwarfs (Garcia-Berro et 
al. 1997a) and. consequently, was not taken into account. 
Second, for the high-density regime, we adopted the equa- 
tion of state described in Segretain et al. (1994) to account 
for all the important contributions for both the liquid and 
the solid phases. In particular, this equation of state, be- 
sides the contribution from partially degenerate electrons, 
accounts for the exchange contribution (Stringfellow et al. 
1990) and the quantum (diffraction) corrections (Hansen 
& Vieillefosse 1975) in the fluid phase, as well as the har- 
monic contribution (Chabrier 1993), the anaharmonic con- 
tribution (Stringfellow et al. 1990), and the electrostatic 
contribution (Stringfellow et al. 1990) in the solid phase. 
For the electrons the exchange effects at finite temperature 
(Kovetz et al. 1972), the correlation contribution (Nozieres 
& Pines 1958), and the electron-ion coupling contribution 
(Yakovlev & Shalibkov 1989) were also taken into account. 
For the low-density regime, we used an updated version of 
the equation of state of Magni & Mazzitelli (1979). Neutrino 
emission rates for pair, photo, plasma, and bremsstrahlung 
processes and for high-density conductive opacities were 
taken from Itoh et al. (1994, 1996a, 1996b). Convection 
was treated in the framework of the mixing length the- 
ory as given by the ML2 parameterization (Tassoul et al. 
1990). Radiative opacities were those of OPAL (Iglesias & 
Rogers 1996), complemented at low temperatures with the 
Alexander & Ferguson (1994) molecular opacities. 

The initial white dwarf models from which we started 
our calculations of the cooling sequences correspond to hot 
white dwarf configurations with a chemical stratification 
appropriate to massive white dwarfs resulting from progen- 
itor stars with solar metallicity that are expected to have 
burned carbon in semidegenerate conditions — see, for in- 
stance Ritossa et al. (1996). According to the calculations 
of these authors, the core is mainly composed of ^^O and 
^"Ne, plus some traces of -'^C, ^'^Na, and ^^Mg and is the 
result of repeated shell flashes that take place during the 
carbon-burning phase in massive intermediate mass stars 
(Garcia-Berro et al. 1997b; Gil-Pons et al. 2003). The core 
chemical profiles of our models are flat throughout the core. 
Such profiles are expected if Rayleigh Taylor instabilities 
act to smooth-out regions with negative molecular weight 
gradients. The outer layer chemical stratification consists 
of a pure hydrogen envelope of 10~^M* (plus a small inner 
tail) overlying a helium-dominated shell of 4 x 10~^M* and, 
below that, a buffer rich in ^^C and ^^O. The amount of 
hydrogen we adopted is an upper limit as imposed by nu- 
clear reactions. All the evolutionary sequences considered 
in this work have the same core composition and shell pro- 
file, which remains fixed during the evolutionary sequences 
and corresponds to that illustrated in Fig. 4 of Corsico et al. 
(2004). The shape of the outer layer chemical profile is given 
by element diffusion at low luminosities. Nevertheless, diffu- 
sion was switched off in the present calculations. Although 
minor changes in the chemical profile are expected because 
of the different masses of progenitor objects, we believe that 
these had a negligible influence on the cooling times of our 
sequences. 

We computed the evolution of our models in a self 
consistent way, using detailed non-gray model atmo- 
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Fig. 1. Evolutionary times as a function of the surface lu- 
minosity for our massive white dwarf models with pure H 
atmospheres. The stages for which neutrino emission, crys- 
tallization and Debye cooling drive the evolution are indi- 
cated. The locii along the curves where the neutrino contri- 
bution to the surface luminosity (100 and 10%), the onset 
of core crystallization (at F = 180) and when the dimen- 
sionless Debye temperature at the core reaches 6-d/T = 20 
are clearly indicated. Note the relatively short ages char- 
acterizing these models at advanced stages of evolution, 
particularly for the most massive ones. 



spheres, which also allow us to derive color indices and mag- 
nitudes of these white dwarfs. Our model atmospheres are 
based on an ideal equation of state. The following species 
have been considered: H, H2, e~, H^, H+, H2+, and H3+. 
All the relevant bound-free, free-free, and scattering pro- 
cesses contributing to opacity were included in our cal- 
culations. At low Teff values, collision-induced absorption 
(CIA) from molecular hydrogen due to collisions with H2 
represents a major source of opacity in the infrared and 
dominates the shape of the emergent spectrum. Collision- 
induced-absorption cross sections are from Borysow et al. 
(1997). Broadband color indices have been calculated using 
the optical BVRI and infrared JHK passbands of Bessell 
(1990) and BesseU & Brett (1988), respectively, with cali- 
bration constants from Bergeron et al. (1997). More specific 
details about the input physics and computational issues 
of the model atmospheres used in this work are described 
at length in Rohrmann (2001). However, for the purpose of 
this paper it is important to realize that in order to compute 
the outer boundary conditions for our evolving models, the 
values of the pressure and the temperature at the bottom 
of the atmosphere are required. In the calculations reported 
here non-gray model atmopheres are fully integrated each 
time the outer boundary conditions must be evaluated as 
evolution proceeds. 

The initial stellar models were derived from an arti- 
ficial evolutionary procedure starting from the full evolu- 
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Fig. 2. Upper panel: The dimensionless Debye tempera- 
ture in terms of the mass fraction for three selected low- 
luminosity 1.28 Mq models. Bottom panel: Internal lumi- 
nosity versus the mass fraction for the same models shown 
in the upper panel. Note that most of the structure of the 
lowest luminosity model is well inside the Debye regime. 



tion of a 0.95 Mq white dwarf model; see Althaus et al. 

(2005) and references cited therein. Despite the correct- 
ness of our artificial procedure to generate starting white 
dwarf configurations, model ages corresponding to the very 
first computed stages of evolution should be taken with 
some care. However, for ages over 10^ yr, the evolutionary 
calculations presented here are independent of the initial 
conditions. We computed the evolution of massive white 
dwarf models with stellar masses of 1.06, 1.10, 1.16, 1.20, 
1.24, and 1.28 Mq. The lower limit of this mass range cor- 
responds to the approximate minimum mass for a white 
dwarf to have an ONe core. The precise value of this mass 
is still not well known. White dwarfs with masses higher 
than ~ 1.1 Mq have progenitors that proceed through car- 
bon burning (Garcia-Berro et al. 1997b), whereas white 
dwarfs with masses lower than ~ 1.0 Mq have progenitors 
that never ignite carbon (Salaris et al. 1997). The evolu- 
tionary stages computed cover the luminosity range from 
log{L/LQi) « 0.5 down to —5.2. 

In addition to the calculation of the evolution of white 
dwarf sequences with pure hydrogen envelopes, we com- 
puted the evolution of the same sequences mentioned above 
but for the case of pure helium envelopes. In this case, we 
assumed a grey atmosphere. Although the employment of 
a gray atmosphere in our helium envelope sequences pre- 
vents a precise quantitative comparison with the case of hy- 
drogen envelope sequences, this still gives us a reasonable 
assessment of the cooling times for these helium envelope 
sequences. 
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Fig. 3. Central temperature versus age and surface lumi- 
nosity as a function of the age (left and right panels, re- 
spectively) for the 1.06, 1.16, and 1.28 Mq sequences. 



3. Evolutionary results 

We begin by examining Fig. [1] which shows the evolution- 
ary times as a function of the surface luminosity for the 
1.06, 1.10, 1.16, 1.20, 1.24, and 1.28 Mq ONe white dwarf 
models with pure H atmospheres. The figure also shows 
(i) the locii along the curves where the neutrino emission 
luminosity becomes 100 and 10% of the surface luminos- 
ity, (ii) the onset of core crystallization (at T — 180), 
and (iii) the dimensionless Debye temperature at the core 
when it reaches 9y}/T= 20, where the Debye temperature is 
Od = 3.48 X 10^(Z/A)p^/^. Several physical processes that 
take place along the evolution leave their signatures in the 
cooling curve. During the first stages of evolution, the in- 
terior temperature is relatively high (between T ~ 1.2 and 
1.5 X 10^ K, depending on the mass of the white dwarf) and 
neutrino emission constitutes a powerful sink of energy that 
considerably affects both the cooling timescales of massive 
white dwarfs and, as shown in Althaus et al. (2005), also the 
mass-radius relation of these stars. Indeed, neutrino emis- 
sion luminosity far exceeds the photon luminosity during 
the hot white dwarf stages. As neutrino emission gradually 
decreases, this accelerated cooling phase arrives at its end 
and the slope of the cooling curve changes, as can be seen 
in Fig. [TJ Note that neutrino losses persist further down 
to luminosities log(L/LQ) ~ —1, but their effects on the 
cooling curve are negligible at those evolutionary stages. 

The release of latent heat upon crystallization is known 
to influence the evolutionary times of white dwarfs as well. 
However, in the case of massive ONe white dwarfs, the crys- 
tallization of the core is a process that occurs at relatively 
high luminosities — log(L/LQ) = —2.15 and —1.30 for 
the 1.06 and 1.28 Mq cooling sequences, respectively (see 
Fig- HI) — and its impact on the evolution of the star there- 
fore remains moderate. It is thus important to note that 
crystallization profoundly influences the cooling behavior 
of only the less massive sequences. 
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Fig. 4. The run of temperature in terms of the outer 
mass fraction for the 1.28 Mq model at luminosities 
log(L/L0) = -3.25, -4.05, -4.50, and -4.96 (from top to 
bottom). Dotted lines represent convectively unstable lay- 
ers. 



The short ages that characterize our white dwarf mod- 
els at very low surface luminosities is a noteworthy feature. 
This is particularly noticeable for our most massive mod- 
els. For instance, the 1.28 Mq cooling sequence reaches a 
luminosity of log(L/L0) = —5.15 in only 4 Gyr. In fact our 
results indicate that the most massive white dwarfs could 
be, depending on the distance, unobservable at ages well 
below 10 Gyr with the current observational facilities. For 
instance, if there are such massive cool white dwarfs at the 
distance of say the Hyades, the apparent magnitude would 
only be my ~ 21.5 or so, and at 1 kpc, mv ^ 28, whereas 
the X-band apparent magnitudes would be 1 or 2 magni- 
tudes brighter, according to our models. At these distances 
such massive cool white dwarfs might be observable with 
present or near future facilities but the observations would 
be very challenging. At the lowest computed luminosities, 
our massive white dwarf models experience the develop- 
ment of the so-called fast Debye cooling. In this sense, we 
would like to mention at this point of the discussion that, 
for the 1.28 Mq model at log(L/LQ) = -4.95, the dimen- 
sionless Debye temperature remains well above 10 for about 
the 95% of the structure (that is, most of the star is within 
the Debye regime), with the consequence that the thermal 
content goes rapidly to zero in this region. For less massive 
CO white dwarfs, this takes place at luminosities of about 
log(L/LQ) = -6 (D'Antona & MazziteUi 1989). This fact 
is reflected in Fig. [21 which displays the run of the Debye 
temperature and the fractional luminosity for some selected 
low-luminosity stages of the 1.28 Mq model sequence. Note 
that for the coolest models, the luminosity output is almost 
negligible in the inner regions of the star and that only in 
the outermost layers are the ratio Oy^/T and the degeneracy 
parameter not very large. Hence, these layers can still be 
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Fig. 5. Same as Fig. [T] but for white dwarf models with 
pure He atmospheres. 



compressioually heated, making the main energy source of 
the star. In fact, the compression of the outermost layers 
provides the bulk of the star luminosity; see Isern et al. 
(1998b) for a detailed description of the cooling of white 
dwarfs. 

For a better understanding of the physical processes oc- 
curring in the interior of our massive white dwarf models 
we show in Fig.[3]thc run of the central temperatures versus 
ages and surface luminosities (left and right panel, respec- 
tively) for the 1.06, 1.16 and 1.28 Mq model sequences. The 
changes in slope of the logTc versus log(L/LQ) diagram at 
both high luminosities and at the low luminosity end are 
noticeable. At high luminosity, it reflects the end of neu- 
trino dominated regime. At low luminosity, the increase in 
the slope of the curve occurs when the hydrogen-rich outer 
convection zone approaches the isothermal degenerate core 
(see Fig. SJ. As a result, the white dwarf has additional 
thermal energy to radiate (D'Antona & Mazzitelli 1989). 
This helps to understand the lengthening of the evolution- 
ary cooling times occurring at log(L/LQ) w —4 in Fig.[TJ 

Non-negligible differences in the cooling of white dwarfs 
arise from the different thicknesses of the H envelopes with 
which white dwarfs proceed during their cooling track. To 
assess such differences, we considered the evolution of mas- 
sive white dwarfs for the extreme situation of pure helium 
envelopes. The resulting evolutionary times are displayed 
in Fig. [5l At advanced stages in the evolution, the central 
temperature of the models is strongly tied to the details of 
the outer layer's chemical stratification. This fact starts to 
manifest itself when the boundary of the degenerate core 
reaches the base of the H envelope at log(L/LQ) ^ —2.5; 
see Tassoul et al. (1990) for details. As a result, white 
dwarfs with helium envelopes (and hence more transpar- 
ent) evolve more rapidly than those white dwarfs with H 
envelopes, as is clear by examining Fig. [S] As expected, the 
helium sequences will reach the Debye cooling conditions 
earlier than their hydrogen counterparts (compare Figs. [1] 




V-K 

Fig. 6. {B — V,V ~ K) color-color diagram for our mas- 
sive white dwarf cooling sequences with masses 1.06, 1.10, 
1.16, 1.20, 1.24, and 1.28 Mq (from top to bottom). Lines 
tracing the locii of equal effective temperatures are labelled 
with the corresponding values. For the 1.06 and 1.28 Mq 
sequences, the ages at Tes = 3600 and 4500 K, respectively, 
are also indicated. Note that all the cooling sequences ex- 
hibit a pronounced turn-off at advanced stages of evolution. 
Below 4500 K, massive white dwarfs become markedly blue 
in this diagram within cooling times much less than 10 Gyr. 



and[5]). Note that the 1.28 Mq white dwarf coohng sequence 
with a helium atmosphere needs only 2.24 Gyr to reach 
log(L/LQ) = —4.70, a factor about 1.5 less than the age 
required by the pure H counterpart. At ages below 5 Gyr, 
most of our cooling sequences with pure He atmospheres 
will have cooled down to below log(L/LQ) « —5.5. 

The fast evolution of our massive white dwarfs at low 
surface luminosities raises the possibility that these white 
dwarfs with pure H atmospheres reach the turn-off in their 
colors and become blue afterwards within relatively short 
ages. To elaborate this point further, we show in Fig. [S] 
the evolution of our pure H white dwarf cooling sequences 
in the {B — V,V — K) color-color diagram. As a result of 
the collision-induced absorption from molecular hydrogen, 
a process that reduces the infrared flux and forces radiation 
to emerge at higher frequencies, very cool white dwarfs are 
expected to become bluer as they age (Hansen 1998). As 
shown in Fig. [6l this effect is also present in our massive 
white dwarf sequences. Indeed, all our sequences exhibit 
a pronounced turn-off in their color indices at advanced 
stages of evolution and become bluer with further evolu- 
tion. In particular, below 4500 K, massive white dwarfs 
become markedly bluer in the {B — V,V — K) two-color 
diagram. The remarkable point is that the turn to the blue 
happens within cooling times much shorter than 10 Gyr. 
If we extrapolate our results, we should thus expect dim 
massive white dwarfs characterized by relatively short ages 
to exhibit blue colors. 
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Table 1. Selected stages for 1.06, 1.16, and 1.28 Mq white dwarf models with H atmospheres. 





Teff (K) 


log(i/Lo) 


Age (10** yr) 


U -V 


V-R 


V-I 


V-K 


B-V 


J-H 


Mv 


1.06 


90000 


0.55 


0.0006 


-1.66 


-0.15 


-0.35 


-1.04 


-0.35 


-0.15 


9.86 


1.06 


85000 


0.44 


0.0010 


-1.65 


-0.14 


-0.35 


-1.03 


-0.34 


-0.15 


9.93 


1.06 


80000 


0.33 


0.0015 


-1.65 


-0.14 


-0.35 


-1.03 


-0.34 


-0.15 


9.99 


1.06 


75000 


0.22 


0.0021 


-1.64 


-0.14 


-0.34 


-1.02 


-0.34 


-0.15 


10.06 


1.06 


70000 


0.09 


0.0030 


-1.63 


-0.14 


-0.34 


-1.02 


-0.34 


-0.15 


10.13 


1.06 


65000 


-0.04 


0.0042 


-1.62 


-0.14 


-0.34 


-1.01 


-0.34 


-0.15 


10.21 


1.06 


60000 


-0.18 


0.0060 


-1.61 


-0.14 


-0.34 


-1.01 


-0.34 


-0.14 


10.30 


1.06 


55000 


-0.34 


0.0091 


-1.60 


-0.14 


-0.33 


-1.00 


-0.33 


-0.14 


10.39 


1.06 


50000 


-0.51 


0.0163 


-1.59 


-0.14 


-0.33 


-0.99 


-0.33 


-0.14 


10.50 


1.06 


45000 


-0.70 


0.0450 


-1.57 


-0.13 


-0.32 


-0.98 


-0.33 


-0.14 


10.62 


1.06 


40000 


-0.91 


0.1655 


-1.54 


-0.13 


-0.32 


-0.96 


-0.32 


-0.14 


10.79 


1.06 


35000 


-1.15 


0.4139 


-1.48 


-0.13 


-0.30 


-0.93 


-0.31 


-0.14 


11.03 


1.06 


30000 


-1.42 


0.8407 


-1.38 


-0.12 


-0.27 


-0.86 


-0.30 
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The molecular hydrogen absorption at low effective tem- 
peratures also affects the evolution of our models in the 
color-magnitude diagram, as shown in Figs. [7] and [3 which 
display the run of the absolute visual magnitude My in 
terms of the V-K and V — I color indices, respectively. 



For the evolutionary stages computed in this work, the turn 
to the blue is noticeable for the V-K color index. Note 
that in this diagram, all our sequences are expected to be- 
come markedly blue for relatively short ages. Specifically, 
our sequences have cooling ages between 3.7 and 6.8 Gyr 



Table 1. Continued. 
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(depending on the stellar mass value) at the turn-off point, 
which occurs at My w 17. For later stages, Debye cooling is 
dominant and evolution indeed proceeds very quickly. Note 
as well that in the {Mv,V — K) diagram our sequences 
remain brighter than My ~ 18 at advanced stages. For 
the lowest luminosities we have driven our cooling tracks, 
our sequences have not yet reached the turn-off point in 
the {Mv,V — /) diagram, as inferred from Fig. [5] For 
the turn to the blue to occur, the evolution should have 
proceeded to lower effective temperatures than those com- 
puted here. Although this would certainly take place in ages 
shorter than 10 Gyr, the expected surface luminosity of the 
models would be extremely low to be detected — below 
log(L/L0) = -5.5. 

In Fig. [HI we illustrate the run of the absolute visual 
magnitude as a function of the effective temperature for 
our massive pure H white dwarf sequences. In addition, we 
draw various isochrones at ages of 0.01, 0.1, 1, 2, and 4 
Gyr. It is clear that the coolest white dwarf models are not 
necessarily the oldest. For instance, in 4 Gyr the 1.28 Mq 
models have cooled down to 4600 K, while the 1.06 Mq 
model sequence remains much hotter (7000 K) at this age. 
Observational data for massive white dwarfs with H-rich 
atmospheres taken from Bergeron et al. (2001) and Liebert 
et al. (2005) are also included in the figure. As compared 
with the ages quoted by the mentioned authors, our results 
suggest younger ages for the white dwarfs in Fig. [SI For 
instance, for WD 1658+441, which, according to our cal- 
culations would have an appreciable fraction of its core in 
a crystallized state, we derive an age of 0.28 Gyr, as com- 
pared with the 0.38 Gyr quoted by Liebert et el. (2005). 
The shorter ages are due partly to the abundant ^°Ne in 



the core of our massive models, which results in a lower 
specific heat per gram. 

Finally, in Fig. [TO] we compare our 1.10 and 1.20 Mq 
sequences with helium atmospheres with those of Garcia- 
Berro et al. (1997a) for ONe cores. Note that our calcu- 
lations predict much younger ages than those derived in 
Garcfa-Berro et al. (1997a). In part, the differences be- 
tween both sets of calculations have their origin in the dif- 
ferent chemical abundance distributions characterizing the 
pertinent models. For instance, our models are character- 
ized by a helium-dominated shell that is four times more 
massive than considered in Garcia-Berro et al. (1997a). It 
should also be taken into account that the preliminary cal- 
culations of Garcia-Berro et al. (1997a) were done using 
a simplified model that obtain the evolution from calcu- 
lating the binding energy very accurately (actually using 
the same equation of state employed here) and coupling it 
with a relationship linking the surface luminosity with the 
core temperature of an otherwise typical CO white dwarf 
of mass 0.6 Mq. 

In closing, we list the main characteristics of our 1.06, 
1.16, and 1.28 A/q white dwarf sequences in Table 1. 
Specifically, we list the effective temperature, the surface 
gravity, the age, the colors, and the absolute visual magni- 
tude. 



4. Conclusions 

We have computed the evolution of massive white dwarf 
stars with oxygen-neon cores for masses ranging from 1.06 
to 1.28 Mq, which covers the expected range of masses for 
which these stars should presumably exist. The evolution- 
ary tracks for both pure H and pure He envelopes were 
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Fig. 7. Absolute visual magnitude My in terms of the color 
index V — K for our massive white dwarf sequences with 
masses 1.06, 1.10, 1.16, 1.20, 1.24 and 1.28 (from top 
to bottom). On each curve, filled symbols represent cooling 
ages, as indicated within the inset. The pronounced turn-off 
at advanced stages of the evolution is apparent. Note that 
for ages well within 10 Gyr, all our sequences are expected 
to become markedly blue in this diagram. 




Fig. 8. Same as Figl7]but for the color index V — I. Unlike 
the V — K, this index has not yet reached the turn-off point 
even for the coolest computed models. 
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computed for surface luminosities spanning the range from 
log(L/L0) « 0.5 down to —5.2. The use of detailed non- 
gray model atmospheres provides us with accurate outer 
boundary conditions for our evolving models at low effec- 
tive temperatures. To our knowledge this is the first at- 
tempt to compute the evolution of massive white dwarfs 
with a realistic equation of state — which includes all the 
non-ideal, corrective terms, and the full temperature de- 
pendence — and reliable chemical profiles for the degener- 
ate interior expected from the previous evolutionary history 
of massive white dwarf progenitors that burned carbon in 
semidegenerate conditions. We have examined the cooling 
ages, colors, and magnitudes of our sequences and find that 
massive white dwarfs are characterized by very rapid evo- 
lution. Indeed, at still observable luminosities, we find that 
the cooling of massive white dwarfs is largely dominated by 
Debye cooling with the result that these white dwarfs could 
be unobservable at ages below the age of the Galactic disk 
with the current observational facilities. At such advanced 
stages, we find our sequences to have reached the turn-off 
in their colors and thus become blue in short times. 

The results presented here will be helpful in interpret- 
ing recent observations of white dwarfs with very high 
surface gravities (Dahn et al. 2004; Madej et al. 2004; 
Nalezyty & Madej 2004), which up to now rely on pre- 
vious evolutionary sequences that were computed assum- 
ing carbon-oxygen cores. Finally, we prepared detailed 
tabulations of ages, colors, and magnitudes for all our 
white dwarf sequences, which arc available at our web site: 
|http : //www ■ f caglp . unlp . edu . ax/ evolgroup/, 
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Fig. 9. Absolute visual magnitude as a function of the effec- 
tive temperature for our massive pure H white dwarf cool- 
ing sequences with masses of (from top to bottom) 1.06, 
1.10, 1.16, 1.20, 1.24, and 1.28 Mq. The thin linos plot 
isochroncs at the labelled ages. Observational data for some 
massive H-rich white dwarfs from Bergeron et al. (2001) 
and Liebert et al. (2005) are included — squares and stars, 
respectively. 
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Fig. 10. Comparison between our 1.10 and 1.20 Mq se- 
quences (helium atmospheres) with those of Garcia-Berro 
et al. (1997a). 
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